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^ Abstract. In this paper we develop and apply methods for the spectral analysis of non-self-adjoint 

P^ tridiagonal infinite and finite random matrices, and for the spectral analysis of analogous deterministic 

r^ matrices which are pseudo-ergodic in the sense of E. B. Davies (Commun. Math. Phys. 216 (2001), 

. 687—704). As a major application to illustrate our methods we focus on the "hopping sign model" 

< l introduced by J. Feinberg and A. Zee (Phys. Rev. E 59 (1999), 6433-6443), in which the main objects of 

" *~| study are random tridiagonal matrices which have zeros on the main diagonal and random itl's as the 

^ other entries. We explore the relationship between spectral sets in the finite and infinite matrix cases, 

r^ and between the semi-infinite and bi-infinite matrix cases, for example showing that the numerical 

I I range and p-norm e-psoudospectra (e > 0, p S [l,oo]) of the random finite matrices converge almost 

surely to their infinite matrix counterparts, and that the finite matrix spectra are contained in the 

T— H infinite matrix spectrum E. We also propose a sequence of inclusion sets for E which we show is 

^ convergent to E, with the nth element of the sequence computable by calculating smallest singular 

^^ values of (large numbers of) n X n matrices. We propose similar convergent approximations for the 

f^^ 2-norm e-pseudospectra of the infinite random matrices, these approximations sandwiching the infinite 

^_^ matrix pseudospectra from above and below. 

1^ Mathematics subject classification (2000): Primary 47B80; Secondary 47A10, 47B36. 

^^ Keywords: random matrix, spectral theory, Jacobi matrix, operators on i^ . 



1 Introduction 



5-H In the last fifteen years there have been many studies of the spectra and pseudospectra of infinite 

random tridiagonal matrices in the non-self-adjoint case, and of the relationship of the spectral 
sets of these infinite matrices to those of corresponding large finite random nx n matrices (see e.g. 
[21, 16, 32, 17, 19, 39, 13, 22, 40, 31, 30] and the references therein). In this paper we contribute to 
this literature, introducing new methods of analysis and computation with emphasis throughout, 
as a major case study, on applying these techniques to understand the "hopping sign model" 
introduced by Feinberg and Zee [17], further studied in Holz, Orland and Zee [22], by ourselves 
previously in [5], and see also [11, 12] and [40, Section 37]. In this model the main object of study 
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is the order n tridiagonal matrix given, for n > 2, by 
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6i 1 



At 







V 



6„-i / 



where 5 = (6i, . . . , 6„_i) € C""^ and each bj = ±1. (For n = 1 we set Af, = (0).) 




Figure 1: A plot of specA^, the set of eigenvalues of Aj^, for a randomly chosen b 6 {itl}"~^, with n = 5000 
and the components bj of b independently and identically distributed, with each bj equal to 1 with probability 1/2. 
Note the symmetry about the real and imaginary axes by Lemma 3.4 below, and that the spectrum is contained in 
the square with corners at ±2 and ±21 by Lemma 3.1 below. 



The objectives we set ourselves in this paper are to understand the behaviour of the spectrum 
and pseudospectrum of the matrix A'^^, the spectrum and pseudospectrum of the corresponding 
semi- infinite and bi- infinite matrices, and the relationship between these spectral sets in the finite 
and infinite cases. Emphasis will be placed on asymptotic behaviour of the spectrum and pseu- 
dospectrum of the finite matrix A„ as n — > oo, and we will be interested particularly in the case 
when the bj are random variables, for example independent and identically distributed (iid), with 
Pr(6j = 1) = 0.5 for each j. (A visualisation of spec A^^ for a realisation of this random matrix with 
n = 5000 is shown in Figure 1; cf. [17].) To be more precise, we will focus on the case when the 
vector b G {±1}"~^ is the first n — 1 terms of an infinite sequence (6i, 62, ■ • ■), with each bj = ±1, 
which is pseudo-ergodic in the sense introduced by Davies [14], which simply means that every 
finite sequence of ±l's appears somewhere in (61, 621 • • ■) a-s a consecutive sequence. If the bj are 
random variables then, for a large class of probability distributions for the bj , in particular if each 
bj is iid with Pr(5j = 1) G (0, 1) for each j, it is clear that the sequence (fei, 62, . . .) is pseudo-ergodic 
almost surely (with probability one). Thus, although pseudo-ergodicity is a purely deterministic 
property, our results assuming pseudo-ergodicity have immediate and significant corollaries for the 
case when Aj, is a random matrix. 



Our interest in studying this problem is in making a contribution to the understanding of the 
relationship between the spectral properties of finite random matrices and corresponding infinite 
random matrices in the difficult non-normal case. (We note that A^ is self-adjoint only in the 
special case that each bj = 1, and it is an easy calculation that A^ is normal, i.e. A^ commutes 
with its transpose, only if 61 = 62 = ■ • • = ^n-i-) For an interesting introduction to the behaviour 
of random matrices in the non- normal case see [40] . Our focus in this paper is on the particular 
matrix A^ and especially its infinite counterparts, but in the course of this investigation we develop 
and apply methods applicable to the study of spectral sets for the much larger classes of infinite 
tridiagonal or banded matrices. 

Our study of the particular matrix A^^, with each bj = ±1, is motivated by interest expressed 
in this class of random matrix in the physics literature [16, 22, 11, 12]. Despite this interest there 
are so far no rigorous mathematical results on the behaviour of the spectrum of A'^^ in the limit as 
n — ^ c». This paper makes steps in this direction. A further motivation for studying the particular 
matrix class A'^^ is that rigorous results are available on the asymptotics of the spectrum and 
resolvent norm for a related class of matrices, offering some hope that progress might be possible 
in this case also. This related class is the case when, rather than the first sub-diagonal consisting 
of random ±l's, the diagonal has random ±l's. Of course, the matrix is then upper-triangular, 
so that many computations become explicit; in particular the spectrum of the finite matrix is just 
{1,-1} and the spectra of the corresponding infinite matrices can be explicitly calculated: see 
[39, 8, 29] for details. We shall see that the situation in the case studied in this paper is, in a 
number of respects, rather richer and the analysis more delicate. At the same time in a number 
of respects our results are more complete: for example, we are able to prove convergence of the 
pseudospectra of A^^ to those of the corresponding infinite matrices, and to do this not just in a 
Hilbert space setting but in p-norm for p G [1, 00]. 

The distinctive flavour of the results we develop in this paper, with their significant emphasis 
on pseudospectra and the relationship between finite random matrices and their infinite matrix 
counterparts, is in large part inspired by the paper by Trefethen, Contedini and Embree [39], by 
Part VI on random matrices in [40] , and by results on convergence of the p-norm pseudospectra 
(1 < p < 00) and numerical ranges of n x n Toeplitz matrices due to Bottcher, Grudsky, and 
Silbermann [1] and Roch [36] , described recently in the monograph of Bottcher and Grudsky [2] . 

Let N denote the set of positive integers and Z the set of integers. Throughout, {±1}^, {±1}^, 
and {±1}™, for m S N, denote the sets of vectors in i°°{Z), £°°{'N) and C™, respectively, whose 
entries bj = ±1. The related infinite-dimensional operators we study include the operators A^, 
for b = (61, 62, . . .) e £°°(N), especially when each bj = ±1. Here A''_^_ acts on the sequence space 
iP{N), for p e [1, 00], by the action 



(A^x), = V(A^).,x„ zeN, (1) 






where {A'^)ij = bi-i6i-i,j + Si+ij and Sij is the usual Kronecker delta. In other words, A'^ acts 
by multiplication by the infinite matrix 
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which has entry {A^)ij in row i, column j, for i,j E N. (For simplicity, we make no distinction 
in our notation between A'^ and its matrix representation.) A main aim of the paper will be to 
compute the spectrum, pseudospcctrum, and numerical range of A'^ in the case when b S {±1}'^ 
is pseudo-ergodic. We shall also study the same properties of the corresponding operator A^ which 



acts on iP{Z), again focusing on the case when b e {±1}^ is pseudo-ergodic. The action of A'' is 



given by the same formula (1) but now with h £ 1°" 
A^ acts by mutiphcation by the bi-infinite matrix 
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(Z) and with N replaced by Z. In other words, 
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where the box marks the matrix entry at (0, 0). Our results will also apply, through the application 
of similarity transforms, to the more general matrices 
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in the case when hj = ±1 and Cj = ±1, and to the corresponding infinite matrices 
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(4) 



1.1 The Main Results 

Let us summarise the main results that we obtain in this paper, first introducing a few key notations 
and definitions. Throughout, where B is a bounded linear operator on (P(S), for some p £ [l,oo], 
with 5 = Z or N, or where i? is a square matrix, we denote by spec-B the spectrum of i3, i.e. 
the set of A G C for which B — XI (/ the identity matrix or operator) is not invertible. (We 
note that the spectra of A'' and A*^ do not depend on p G [l,oo], from general results on band 
operators (e.g. [26]); of course, when B is a matrix, the spectrum is just the set of eigenvalues 
of B.) Throughout, ||a;||p, for p G [l,c»], will be our notation for the standard p-norm of x, for 
X G P'iS), with 5 = Z or N, or a; G C™, for some m G N. Where B is an operator or matrix, \\B\\p 
will denote the norm of B induced by the vector norm || • ||p, i.e. \\B\\p := supiun ^^ ||i?x||p. With 
this notation, following e.g. [40], for p G [1, 00] and e > we define the P e-pseudospectrum of B, 
specfi?, by 

specks :^ specs U {A gC : \\{B - \iy^\\p > e-^}. 

When B is a bounded linear operator on ^{S), for some p G [l,oo] and S* = Z or N, in general 
the spectrum of B is larger than the set of eigenvalues of B. We let speCp^jj^^ B denote the set of 
eigenvalues of B considered as an operator on P'{S), i.e. 



spec^oint B 



{XeC:Bx = Xx, for some x G F(S') with x ^ 0}. 



A key result we obtain on the spectra of our infinite matrices, in large part tlirougli limit operator 
arguments described in Section 2, is the following (cf. [f4]): if b,c,d £ {±1}'^, b,c,d £ {±1}^, and 
6, &, cd, and cd are all pseudo-ergodic, then 

spec A''^ = spec A"/^ = spec A^ ^ spec A'^'^ = Y. := (J specA'^ = (J spec[^i„tA^ (5) 

ee{±l}S ee{±l}2 

One surprising aspect of this formula is that the semi-infinite and bi-infinite matrices share the 
same spectrum, in contrast to many of the cases discussed in [39] , this connected to the symmetries 
that we explore in Section 3. 

We do not know a simple test for membership of the set S given by this characterisation (though 
see Figures 2 and 3 below for plots of known subsets of E, and see Section 4.3 for an algorithm 
for computing approximations to E). But this result implies that spec^'' C E for every b G {±1}^ 
which gives the possibility of determining subsets of E by computing spec A'' for particular choices 
of b. In particular, as recalled in Section 2, when b is n-periodic for some n G N, i.e. 6j+„ = bj 
for j S Z, specyl'' can be computed by calculating eigenvalues of an order n matrix (a periodised 
version of A^). We compute 7r„ C E, for n = 5, 10, ..., 30 in Section 2, where 7r„ denotes the union 
of spec Af, over all n-periodic b € {±1}^. We speculate at the end of the paper that 

7i"oo := U 7r„ (6) 

is dense in E, and it has been shown recently in [9] that certainly tToo is dense in the unit disc 
D = {z : |z| < 1}, which implies that D C E, as established slightly earlier directly from (5) in [5]. 
(Throughout, S denotes the closure of a set S* C C: for an element z G C, z denotes the complex 
conjugate.) 

To obtain a first upper bound on E we compute the ^^-numerical range, Vl^(A''), of A^ when 
b is pseudo-ergodic. We show that, if b,c,dG {±1}'^, b,c,d E {±1}^, and b, b, cd, and cd are all 
pseudo-ergodic, then 

WiA\^) = WiA"/) = WiA'') = W{A^^"^) = A:={z = a + ib:a,beR,\a\ + \b\ < 2}. 
Since the spectrum is necessarily contained in the closure of the numerical range, this implies that 

D C E c A. 

We point out that the numerical range of A'^^ converges to that of A^, in particular that VT (A^) /^ A 
as n — 7> oo, if 6 is pseudo-ergodic. (Here and throughout, for T„ C C and T C C, the notation 
Tn /• T means that T„ C T for each n and that dist (T,T„) -^ as n ^- oo, with dist(T, r„) the 
Hausdorff distance defined in (16) below.) 

The largest part of the paper (Section 4) is an investigation of the relationship between the finite 
and infinite matrix cases with respect to behaviour of spectra and pseudospectra. The spectral 
case is harder: our main result is to show that the spectra of the finite matrices are subsets of the 
infinite matrix spectra, precisely that, for every n and every c € {±1}"^^, 

specA^ C 7r2„+2 C E, 

so that o-„ := Uce{±i}"-i specA'^ C 7r2„+2 C E and 

CToo := y Cr„ C TToo C E. (7) 

nSN 

We suspect that specA^ /^ E = specA\_ as n — > oo, if 6 e {±1}^ is pseudo-ergodic, and the 
numerical results in Figures 1 and 2, and other similar computations, are suggestive of a conjecture 
that specA'^ /• Wao, which set, as mentioned already, we speculate is dense in E. 



We can prove neither of these last two conjectures about spectral asymptotics. On the other 
hand, our theoretical results for the pseudospectrum are fairly complete. We show first in Theorem 
3.6 a pseudospectral version of (5), that, if 6, c, d G {±1}^, b,c,d G {±1}^, and b, b, cd, and cd are 
all pseudo-ergodic, then, for p £ [1, oo] and e > 0, 

spec^A'; = apecPA"/ = spec^A^ = spec^yl'^''^^ = SJ? :- |J spec^A^ 

ee{±iF 

We then show that the pseudospectra of the large finite matrices are contained in and are well- 
approximated by the pseudospectra of the infinite matrices, and that this works for p-norm pseu- 
dospectra for the full range p G [l,oo]. Precisely, for p E [l,oo] and e > 0, we show that, if b is 
pseudo-ergodic, then 

spec^^:; ^ EP (8) 

as n — >■ cx). 

This last result, linking the pseudospectra of A^ with those of its finite sections A^, is a 
somewhat unexpectedly satisfactory result. Even in the case in which the theory of the finite 
section method is arguably simplest and most well-understood, namely the case of the Toeplitz 
operator (a semi-infinite Toeplitz matrix) , the limit as n — >■ oo of the e-pseudospectra of the n x n 
finite section Toeplitz matrices has been calculated only relatively recently, and only for p G (1, oo) 
[1, 2]. Moreover, except for the special case p = 2, this limit is not, in general, just the i^ e- 
pseudospectrum of the Toeplitz operator, but rather the union of the l^ and i"^ e-pseudospectra, 
with p~^ + q~^ = 1. (A component of the explanation of (8) is that we show in Lemma 3.2 that 
Ef = S9 forp-i +q-^ = 1.) 

Equation (8) leads to characterisations of the spectrum S which, in principle, can be used for 
numerical approximation. Since ne>o ^e ^ ^' ^^ holds that 

E = lim EP = lim lim spec^^^, (9) 

e— fO e— ^0 n— foo 

for every p e [1, oo] and pseudo-ergodic b. However, the formula (9) is not guaranteed to give useful 
results for any fixed e and n as the convergence as n — > oo may be arbitrarily slow, as discussed 
in Section 4.3. In that section we develop alternative, much more useful, convergent sequences 
of computable, upper and lower bounds for E^ and a convergent sequence of computable upper 
bounds for E. We show firstly that 



al 



U spec^A^ c e2 c (jI,+,„ = U spec^+.^A^ 
cG{±i}"-i ce{±i}"-i 



giving explicit expressions for the e„ which satisfy that e„ — 0{n''^) as n — )■ oo, and showing that 
(j^ £ /* E^ and tr^ ^^^ \ ^? as n — >■ oo. (The notation Tn \ T means that T C r„ for each n 
and that dist(r, r„) — > as n — >■ oo.) Then, taking the intersection over all e, we deduce that 



an = IJ spec A^ C E C a^, 
ce{±i}"-i 

and prove that 

(7^ ^ \, E as n — ?► oo. 



In a substantial series of numerical calculations, we compute these convergent upper bounds a'^ ^ 
for the spectrum E in Section 4.3, and through these calculations demonstrate that E is a strict 
subset of A. 

All these results have implications for the behaviour of the spectral sets of A'^, A'^, A'^, A^''^, 
Aj^"^, and A^'^, when the entries bj — ±1 and Cj — ±1 are random, and we make explicit these 



implications in a final Theorem 5.1, in the same section summarising succintly what we have 
established about the spectral sets S and Ef (Theorem 5.2), and outlining a number of open 
problems. 

In the course of this investigation, focused on a particular operator and matrix class, we develop 
results for the larger classes of tridiagonal or banded finite and infinite matrices. In particular, 
Theorem 4.4 shows that, for p S [1, oo], £ > 0, the £p e-pseudospectrum of a general, semi-infinite 
tridiagonal matrix is contained, for e' > e, in the l^ e'-pseudospectrum of its n x n finite section if 
n is sufficiently large. It also shows corresponding results relating the pseudospectra of a general 
bi-infinite matrix to that of its finite sections. In Section 2 we employ recent work [7, 8] on limit 
operator methods for the study of spectral sets for very general classes of infinite matrices. We 
make explicit in Theorems 2.1 and 2.9 the implications of this work for the essential spectrum, 
spectrum, and pseudospectra of bi-infinite and semi-infinite banded matrices with numerical (as 
opposed to operator-valued) entries. In Section 4.3 we make the first substantive application of 
a new method which generates sequences of inclusion sets for the spectra and pseudospectra of 
a tridiagonal operator, demonstrating, through this application, that these sequences of inclusion 
sets can in fact converge to the spectral sets that they enclose. 

We mention that first versions of a number of the results in this paper are contained in the PhD 
thesis of the second author [10], and that a number of the results were announced (without proofs) 
in [5]. 

1.2 Pseudospectra and the Numerical Range 

We shall need throughout the paper a number of properties of the e-pseudospectra of a bounded 
linear operator _B on a Banach space X, and of the pseudospectra of its adjoint operator B* 
on the dual space X* (dual in the sense e.g. of [24], so that X* is the set of bounded anti-linear 
functionals, and the spectrum of -B* is the complex conjugate of the spectrum of -B) . We summarise 
these properties in this section, pointing out how the theory of pseudospectra in the Banach space 
setting has recently been significantly clarified by work of Shargorodsky [37]. The properties we 
shall need include the equivalent definitions encapsulated in the following theorem: 

Theorem 1.1 The e-pseudospectrum of a bounded linear operator B on a Banach space X is 
defined, for e > 0, by any one of the following equivalent definitions: 

(i) specks = spec BU {A eC : ||(S- A/)-i|| > e-^}; 

(ii) speCgi? = speci?U{A e C : v{B — \T) < e}, where v{C) is the lower norm of a bounded linear 
operator C, defined by i'{C) := inf||j.||^i lie's;]!,- 

(iii) specks is the union of spec B and the set spec jj^j ^ B of e-pseudoeigenvalues of B, where A 
is an e-pseudoeigenvalue if there exists x Cz X with \\x\\ ~ 1 and \\(B — A/)a;|| < e; 

(iv) speCgi? is the union o/ speCpoj^t ^ i? and the complex conjugate o/speCpoj„t ^ i?*; 

(v) specks = U||£;||<^ spec (S + E), the union taken over all bounded linear operators E with 

\\E\\<e. 

For a proof of the equivalence of (i)-(v), and a useful short introduction to the pseudospectra 
of linear operators on Banach spaces, see [40, Section 4]. We will use the equivalence of (i)-(iv) 
throughout. The equivalence of the other definitions with (v), and the connection this makes with 
spectra of perturbed operators, is a significant motivation for the practical interest in pseudospec- 
tra. It is clear from the above definition that spec _B is an open set for e > 0. An elementary but 



important property of the lower norm is that 

HA)~v{B)\<\\A-Bl (10) 

for any bounded hnear operators A and B on X. 

In the case when, for some iV e N, X = C^ and i? is an A^ x A^ matrix, (i)-(v) are equivalent 
additionally to spec^i? = {A e C : v{B — XI) < e} = speCpQj„t ^ i?. If || • || = || • ||2, then, for every 
N X N matrix A, lyiA) = s„iin(v4), the smallest singular value of A. Thus these definitions are 
additionally equivalent to [40] 

specks = {A e C : Snii„(B - XT) < e}. (11) 

Note that (10) implies that 

\s^in{B~XI)~Sn,in{B-flI)\<\X-fl\, X,fieC. (12) 

It is equation (11) that we use for the numerical computations of pseudospectra in Section 4.3. 

An alternative definition of the pseudospectrum is to replace the strict inequality > in (i) by 
>, so that the £-pseudospectrum is defined to be 

Spec^B = specSU{AeC: \\{B - XI)-^\\ > e"^}. 

This has the attraction that Spec^B, like spec B, is a compact set for e > 0. An interesting question 
is whether spec^-B — Spec^B, which hinges on the question of whether or not it is possible for the 
norm of the resolvent of S, || (B — XI)^^\\, to take a finite constant value on a open set G C C. Let 
us say that the Banach space X has the strong maximum property if, for every open set G C C, 
every bounded linear operator i? on A", and every M > 0, it holds that 

(||(B - XI)-^\\ < Af, VA e G) ^ {\\{B - A/)-i|| < M, VA e G). 

If X has the strong maximum property, then no bounded linear operator on X can have a resolvent 
norm with a constant finite value on an open subset of C, and it is easy to see that spec^i? = Spec^B. 
Recently, Shargorodsky [37] has shown, by constructing explicit counterexamples, that not every 
Banach space has the strong maximum property. But the following theorem from [37], which 
extends earlier work of [18], makes clear that the Banach spaces of relevance to this paper do have 
this property. 

Theorem 1.2 Suppose that X is a Banach space which is either finite- dimensional or is such 
that either X or X* is complex uniformly convex (as defined e.g. in [37]). Then X has the strong 
maximum property. In particular, X has the strong maximum property if X is a Hilbert space, or 
ifX = iP{S),forS^N orZ and p e [l,oo]. 

It is clear from (v) and standard operator perturbation arguments (see [40] for details) that, for 
< e < e', speci? C spec^i? C spec^,B, and that 

eD + spec B c spec^B. (13) 

In fact [40] eD + specB = spec^i? if A" is a Hilbert space and B is normal, i.e. BB* = B*B. 
Further [40] 

speci? — ll specks. (14) 

Generalising (13), it holds that [40] 

(5D + spec^B c spec^+^B, for e,6 >0. (15) 



For 5,rcC, let 

dist(S', T) := max(sup{dist(z, S) : z e T}, sup{dist(z, T) : z e S}). (16) 

(This notion of distance, when apphed to compact subsets of C, is an instance of the Hausdorff 
distance between compact subsets of a metric space.) Given a sequence T„ C C and T C C, let 
us write r„ -;► T if dist(T„, T) ^ asn-^ oo. Additionally, let us write r„ ^ T if T^ -> T and 
Tn C T for each n, and write Tn\T ii Tn ~^ T and T C T„ for each n. It is an easy calculation 
to show that 

speCgi? \ specB as e — )• 0^. (17) 

Similarly, it holds for e > that spec^,B \ Spec^i?, as e' — )■ e+, and spec^,B /• apec^B, as 
e' — > e~. Thus, in the case where X has the strong maximum property so that spec^i? = Spec^JS, 
it holds for e > that 

spec^/i? \ speCgi?, as e' — ^ e^, and spec^/i? /* spec^B, ase'— >-e~, (18) 

so that speeds depends continuously on e. 

The spectrum and e-pseudospectra are connected to the numerical range. In the case that X 
is a Hilbert space with inner product (•, •), and where i? is a bounded linear operator on X, the 
numerical range or field of values of B, denoted W{B), is the set 

W{B) := {{Bx,x) -.xeX, \\x\\ = 1}. 

It is well known that this numerical range is a convex set and that specB C W(B), in fact 
specB C W{B) if X is finite-dimensional. The relationship with the e-pseudospectra is that, 
similarly, spec^i? C W{B) + eD, for e > [40, Section 17]. Let F be a closed subspace of X, 
P : X ^ Y orthogonal projection onto Y, and let By '■= PB\y- Then 

W{By) = {{Byx,x) -.xeY, \\x\\ = 1} = {iBx,x) -.xeY, \\x\\ = 1} C W{B). (19) 

This observation is one component in the following result [20, Theorem 3.52]: 

Theorem 1.3 Suppose that X is a Hilbert space and that (P„)„gN is a sequence of orthogonal 
projection operators on X that converges strongly to the identity operator (PnX -^ x as n — > cxd, 
for every x S X). Then, for every bounded linear operator B on X , where _B„ := P„_B|jf„ with 
Xn — Pn{X), it holds that 

W{B,,) /" W{B) as n -> oo. 

2 Results by Limit Operator Arguments 

Let us start this section by establishing a few additional notations and definitions. Throughout 
the remainder of the paper, if _B is a bounded linear operator on a Banach space X we will say 
that B is Fredholm if B(X), the range of i?, is closed and if, additionally, a(i?) := dim(keri3), 
the dimension of the null-space of i?, and /3(-B) :— dim(X/i?(X)), the co-dimension of the range 
of B, are both finite, in which case we define the index of B by indS := a(B) — P{B). We will let 
speC(,gg B denote the essential spectrum of B, i.e. the set of A € C for which B — XI is not Fredholm. 
Let Mb be the bounded linear operator which operates on the standard sequence space £p(Z), for 
p e [1, oo], by multiplication by 6 e £°°{Z). ExpHcitly, for y e F(Z), 

{Mby)j=bjyj, jeZ. 

Moreover, for A: e Z let Vk denote the shift operator defined by 

{Vky)j = Vj-k, j e Z, 



and note that VjMi, = My bVj, for j E Z, b E £°°{Z). In terms of these notations, the operators 
A^ and A''''^, corresponding to the infinite matrices (2) and (4), can be written as 

A'' = ViMi, + V^i and A''^" = ViMf, + McV^i. (20) 

We wiU use these notations for b,c E i°°{Z), but especially for b,c E {±1}^. 

One major tool for computing the spectrum of the infinite matrices A'' and A'^''^, with b,c E 
£°°{Z), is the method of so-called limit operators [8, 26, 35]. In this method a bi- infinite matrix B 
is studied in terms of a family of infinite matrices that represents the behaviour of B at infinity. 
More precisely, let A be a banded matrix A ~ (ay)ijgz, with supj^ \aij\ < oo, so that the operator 
induced by ^ is a bounded operator on iP{Z), for all p E [1, oo]. We say that the operator induced 
by the matrix B = (6ij)ijgz is a limit operator of the operator induced by A if, for a sequence 
hi,h2, ■■■ of integers with \hk\ -^ oo, it holds that 

o-i+hkj+ht -^ hj as k ^ CO, 

for all i,i E Z. The set of all limit operators of A is denoted by cr°P(A). In some instances it is 
useful to think of cr°P(A) as the union of two subsets, as (t°p(A) = (7°^ {A) U a°^{A), where (t^{A) 
denotes the subset of those limit operators associated with sequences h with hk — > ±oo. It is 
an easy consequence of the Bolzano- Weierstrass theorem and a diagonal argument that each of 
a°^{A) is non-empty, and it is clear that if i? = {bij) is a limit operator of A then sv^Y>i-j=k \^ij\ ^ 
supj_^^^ \aij\, for every k E Z. In particular, ii A = A'^'^ for some b,c E {±1}^ and _B is a limit 
operator of A, then B = A^''^ for some b,cE {±1}^. 

The following theorem, which applies in particular to A^ and to A^-'^, connects the essential 
spectrum with the set of limit operators. This result is a particular case of much more general 
results from [7], [8, Theorem 6.28, Corollary 6.49], which extend a main theorem on limit operators 
going back to [25, 34]. Note that the spectrum, as an operator on £p(Z), of an infinite banded 
matrix A = {aij)ij^z, with supj^ \aij\ < oo, does not depend on p G [l,oo], and the same is true 
for the essential spectrum: moreover, if A ^ spec^^gg A, then ind {A — XI) is also independent of p 
(see [28] or [8, Corollary 6.49]). 

Theorem 2.1 Let A be a handed matrix A = {aij)i,j,=z, with sup^ \aij\ < oo. Then 

spec^gg A^ y spec B = [j speCj^i^j B (21) 

Becr°i'{A) BecT°P{A) 

and specfi? C specl*^, for all e > 0, p E [l;Oo], and B E a°'^{A). In particular, if A E (t°p{A), in 
which case we say that A is self-similar, then 

spec A = spec^gg A — M speci?= M spec'^^l^^B and spec^^ = M spec^B, 

BeCT°P(A) Becr'i'iA) Secr°P(A) 

for £ > and p E [I, oo] . 

Recall that specj^jj^j B is the set of eigenvalues of B in £°°{Z), so that A E specj^j^j i? iff Aa; = Bx 
has a non-trivial bounded solution. 

One case where A^'"^ is self-similar is where (6, c) is periodic with some period n G N, i.e. 

bj+n = bj, Cj+n = Cj, j E Z. (22) 

In this case the above theorem applied to A''''^ reduces to spec A''"'^ = speCggg A^''^ = spec^j^^ A^''^, 
and in fact it is well-known further, e.g. [15], that if A G spec^''''^ then Ax = A''''^x has a solution 
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which is not only bounded but also quasi-periodic, i.e. for some a G C with |a| = 1, Xk+n = o.Xk-, 
k £ Z. It is easy to see that this implies that 

spec A^'^ = U spec {A^^^ + B'^%) , (23) 

|a| = l 

where A'^'^ is given by (3) (with A^'^ := (0)) and -B^'^ is the n x n matrix whose entry in row i, 
column j is (5i.„^j_iac„ + (5i_i(5j,„a^^fe„, where Sij is the Kronecker delta. We will abbreviate -B.^'^ 
as B^ ^ in the case that c = (1, ..., 1). 

An important case where A'' is self-similar is where A'' is pseudo-ergodic in the sense of Davies 
[14]. The following is a specialisation of the definition from [14]. 

Definition 2.2 Call b € {±1}^ and the operator A^ pseudo-ergodic if, for every N E N and every 
w e {±1}^, there exists J G Z such that bn+j = Wn, for n = 1, ...,7V. 

We see from this definition that A^ is pseudo-ergodic if and only if every finite sequence of 
±l's appears somewhere in the bi-infinite sequence b. The significance of this definition is that, 
for many cases where the entries 6„ are random variables, the sequence b is pseudo-ergodic with 
probability one. In particular, the following lemma follows easily from the Second Borel Cantelli 
Lemma (e.g. [3, Theorem 8.16]), the argument sometimes called the 'Infinite Monkey Theorem'. 

Lemma 2.3 // the matrix entries 6„, for n S Z, are iid random variables taking the values ±1 
with Pr(6„ = 1) € (0, 1), then A^ is pseudo-ergodic with probability one. 

The link to limit operators is provided by the following lemma (see [14, Lemma 6], [26, Corollary 
3.70] or [8, Theorem 7.6]): 

Lemma 2.4 For b G {±1}^, A'' is pseudo-ergodic if and only ifa°P{A'') = {A'' : c e {±1}^}. 

Combining this lemma with Theorem 2.1 gives the following characterisation of the spectrum 
and pseudospectrum of A^ in the case when b is pseudo-ergodic: 

Theorem 2.5 If b E {±1}^ and A'' is pseudo-ergodic, then 

specA'' = speCe3s^''= |J spec^= = S:= |J spec^,i„t ^" (24) 

ce{±i}z ce{±i}2 

and 

specP A'' == EP := |J specP^^ (25) 

ce{±iF 

for £ > and p £ [I, oo] . 

Limit operator ideas, the "Infinite Monkey" argument and the validity of the first two "=" 
signs in (24) are not new in the spectral theory of random matrices (see e.g. [4, 13, 14, 19, 33]). 
Equation (25) is previously shown, for a general class of pseudo-ergodic operators for the case p = 2 
in [14]. What is more recent is the third "=" sign in the first of equations (24) and the extensions 
to p e [l,oo], these shown in [7] and [8, Theorem 6.28, 7.6]. 

Note that the above theorem shows that the spectrum of A'^ is the same set S for every pseudo- 
ergodic b € {±1}^, and that specA'^ C S for every c € {±1}^, and that similar statements hold 
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for the pseudospectrum specfA^. In particular, specA'^ C S if c G n„, for some n £ N, where 
n„ := {c e {±1}^ : c is n-periodic}. Thus 

7r„ := y spec^^ = |J spec^„i,,t ^^ C E, (26) 

cen„ c6n„ 

for every n E N: this is informative as 7r„ can be computed exphcitly by (23) as the union of 
eigenvalues of n x n matrices. The following lemma carries out this computation for n = 1, 2, 3. 

Lemma 2.6 If belli with bo = I, then spec A'' = [-2,2] and spec^"'' = i[-2,2]. //6 G Ha \ Hi 
then specA'' = T2 '■= {xzL'ix : —1 < a; < 1}. If b Q H^, b^ = bi = 1, and 62 = ^1; then 

spec A^ = T3 := i[-l, 1] U {a; + iy : -1/2 <y < 1/2, a;^ = 1 + 3y^} 

while spec A~'' = iT^ . Thus 

^i = [-2,2]Ui[-2,2], 7r2-7riUr2, TTg = ^1 U T3 U irg. 

Note that max;^g^J |A| — 2 while max>,g,-. |A| = \/2, for j — 2,3. For j = 2 this maximum is 
achieved at ±1 ± i, while for J = 3 this maximum is achieved at ±v7/2 ± i/2. 

Proof, lib eUi with 60 = /3 = ±1 then, from (23), specA'' = Ui^i^ispecSj „ = {e'^ + e^'^/3 : 
ee R}. So spec^'' = [-2,2] if ;3 = 1 and spec A^ == i[-2, 2] if /3 = -l, and tti ^ [-2, 2] U i[-2, 2]. 

If fe e Ha \ Hi then, from (23), where /3 = 61 = ±1, 

spec A^ = y spec ( ^^^g ^ ~ ^^ ' ^ j = {A e C : A^ = 2isin6i, 6* e M}. 

Thus spec A = T2 and 7r2 = vri U r2 . 
If & e Ha, bo = bi = 1, and &2 = 

specv4^= y spec I 1 1 ] = {A e C : A^ - A = -2i 

6»eK y e'^ 

Writing A = a; + iy, we see that A'' — A = — 2isin0, for some 6* G M, iff 

x(a;2 - 3j/^ - 1) = and 3x^y - y^ - y e [-2,2]. 

But this implies that either x = and y^ + y e [-2, 2], or x"^ = 3y^ + 1 and 8y^ + 2y e [-2, 2], and 
it follows that spec^*" — T3. That spccA^^ = it can be shown similarly, or follows from Lemma 
3.4 below. Since c e lis iff c = ±Vjb for j = 0, 1 or 2, it follows that vrs = tti U T3 U ir^. ■ 

In Figure 2 we plot 7r„ for n — 5, 10, ...,30, with 7r„ computed numerically in Matlab using the 
characterisation (23) (see [5] for small plots of 7r„ for n = 1,2, ...,30). For each n the set 7r„, by 
the characterisation (23), consists of A; < n2" analytic arcs, and 7r„ C S. The visual impression 
that might be taken from this sequence of plots is that 7r„ "fills out" a large part of the square 
A := {x + iy : x,2/ G M, [a;] + \y\ < 2} as n — >■ 00. But of course tToo := ^netiT^n is a countable union 
of analytic arcs, so that tToo has (two-dimensional) Lebesgue measure 0. Thus almost every point 
in A is not in tt^o and so is not one of the points in the plots in Figure 2. Thus these figures provide 
no evidence that the Lebesgue measure of E is any larger than zero. And indeed it was conjectured 
in [22] that E has fractal dimension in the range (1,2) (and so Lebesgue measure zero). That this 
is not the case was shown in [5] by an application of Theorem 2.5, specifically by constructing a 
sequence b G {±1}^ for which spec^^jj^j A^ D D, the open unit disc. Of course, this implies by 
Theorem 2.5 the following result. 
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Figure 2: Our figure shows the sets it„, as defined in (26), for n = 5, 10, ..., 30, computed using the characterisation 
(23), which is made explicit for n = 1, 2 and 3 in Lemma 2.6. In particular, tti = [—2, 2] U i[— 2, 2] and, for each n, 
TTi C TTn and, by Lemma 3.1, it„ C A = {a; + iy : i, j/ S K, |a;| + |y| < 2}. 



Theorem 2.7 [5, Proposition 2.1] D C S. 



Recently [9], an alternative proof of this theorem has been obtained, through a construction that 
shows that tt^ is dense in D. It is an open (and interesting) question as to whether t:^ is dense in 
E. An interesting, related, case where the union of the spectra of all periodic operators is shown 
to be dense in the spectrum of the pseudo-ergodic case is studied in [29], but there are other 
pseudo-ergodic bi-infinite tridiagonal examples where this is not true. 

The above results concern bi-infinite matrices, but similar results apply to the semi-infinite 
matrices A5J_ and A_^ . We say that the operator induced by the bi-infinite matrix B = (fey)i.jgN 
is a limit operator of the operator induced by the banded semi-infinite matrix A+ — {aij)ij^fi if, 
for a sequence hi,h2, ■■■ of integers with hk — > +oo, it holds that 



O'i+hk.j+hk 



k — ?► oo. 



for all i,j e Z. The set of all limit operators of A+ is denoted by a°'^{A+). An equivalent 
characterisation is that a°^{A^) = (7°^{A^), where, for any semi- infinite matrix A^, A^ is the 



bi-infinite matrix defined by A-|- = {aij)ij^z, where di 



,-, i,j € N, di 



0, otherwise. The 



following version of Theorem 2.1 holds in the semi-infinite case. In its results on the pseudospectrum 
this theorem appears to be new and may be of independent interest. The arguments in this theorem 
and in later sections depend on the following lemma which, in its results for the pseudospectrum, 
generalises [40, Theorem 2.4(iii)] from the finite-dimensional Hilbert space case to an infinite- 
dimensional Banach space setting, and so may also be of independent interest. 
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Lemma 2.8 Suppose that X is a Banach space which can be written as the direct sum of two closed 
subspaces as X ^ Xi © X2, by which we mean that each x Cz X can be written in a unique way 
as X — xi + X2 with xi S Xi and X2 G X2, and that there exists a continuous projection operator 
Pi : X ^ Xi (in which case P2 = I — Pi is a projection operator onto X2). Suppose also that A is 
a bounded linear operator on X which has Xi and X2 as invariant subspaces, and let Aj denote A 
restricted to Xj, for j = 1,2. Then spec^ — spec^i Uspec A2, speCggg ^4 = spec^^gg ^1 UspeCggg A2, 
and spec^Aj C spec^A, for e > 0, and j — 1,2. If, for some p G [1, 00], it holds for every xi € Xi 
and X2 G X2 that \\xi + X2\\ = ||(||a::i||. ||a;2||)||p, then also spec^^ = spec^Ai U spec^A2, for e > 0. 

Proof. The identities spec A — spec Ai U spec A2 and spec^^gg A — speC(,gg Ai U speC(,gg A2 are 
standard, see e.g. [23, 15]. By Theorem 1.1, spec^i? = speci? U {A G C : i^{B — XI) < e}. 
Since i^iAj — XI) > v{A — XI), for all A G C and j = 1,2, it follows that spec^Aj C spec^A, 
for e > 0, and j — 1,2. II, for some p G [l,oo], it holds for every xi G Xi and X2 G X2 that 
lla^i + X2\\ — IKIIxill, |ja;2||)||p, then, for every A G C, where B := A — XI and Bj := Aj — XI, for 
j = 1, 2, it holds for xi G Xi and X2 G X2 that 

\\B{xi+X2)\\ = \\{\\Bxi\\,\\Bx2\\)\l>\\{HBi)\\xi\\,iy{B2)\\x2\\)\\^ 

so that 

,„, ., \\Bixi+X2)n .^ \\HBi)\\xi\\,iyiB2)\\x2\\)\\ 

i/(B) — mi — r — > nit 



a;iGXi,a;2eX2 ||2:i+a;2|| nGXi ,i;2eX2 || (ll^^i || , ||a;2 1|) || 

But it is an easy calculation that this last infimum has the value min(i^(_Bi), v{B2)). Thus spec^A C 
spec^Ai U spec^A2. ■ 

Theorem 2.9 Let A^ be a semi-infinite banded matrix A^ — {aij)ij^jq, with sup^ \aij\ < 00. 
Then 

speCeggA+= IJ specs = |J spec^^^^^B. (27) 

Further, specf B C specf A+, for all e > 0, p £ [1, 00], and B G ct°p(A+). 

Proof. Given e > and p G [l,oo], choose A > ||A+||p + 2e and apply Theorem 2.1 to the 
bi-infinite matrix A = A^ + B, where B = {bij)ij^i is defined by 6y := X, ii i = j < 0, 6y := 0, 
otherwise. Since a°P{A) = a°^{A) U a°J'{A) = a°P{A+) U {XI}, we see, applying Lemma 2.8, that 

{A} U spec^gg A+ = spec,gg A = {A} U y specB = {A}U |J spec^„i„t B. (28) 

BeCT<'p(A+) _BeCT°p(A+) 

Since A > ||j4+||p > ||-B||p is not in spec^^gg yl+ or in speci?, for B G (t°p(j4+), equation (27) 
follows. Similarly, applying Lemma 2.8, specf A = specf (A/) U specf yl_|- = (A + eD) U spec^A^. It 
follows from Theorem 2.1 that, for B G a°P{A+), specf B C specfA = (A + eD) U specf A+. Since 
specf B c (||B|jp + £)D and A > \\A+\\p + 2e> \\B\\p + 2e, this implies that spec^B c specfA+. ■ 

One consequence of this result and Theorem 2.1 is the following lemma. 

Lemma 2.10 For b,c£ {±1}^, G spec^gg A'''" and G spec^gg A^ ■ 

Proof. It is easy to see that G specS^jj^j A'^''^ for every 6, c G {±1}^, and the result then follows 
from equations (21) and (27). ■ 

We extend the definition of pseudo-ergodic in Definition 2.2 to the semi- infinite case in the 
obvious way, replacing Z by N and A^ by A*^, so that b G {±1}^ is pseudo-ergodic iff the sequence 
b contains every finite pattern of ±l's. Then Lemma 2.3 holds with Z replaced by N and A'' 
replaced by A^, and Lemma 2.4 holds with A^ replaced by A^. 
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3 The Numerical Range and Symmetry Arguments 

Let us first introduce some properties of and notation related to adjoint operators. Given a banded 
bi-infinite matrix A = {aij)ij^z, witli supj- \aij\ < cx), A* will denote the matrix A* = {aji)i.j(zz- 
For 1 <p < oo, where q e (1, oo] satisfies p^^ + q^^ = 1, and identifying ^"^(2) with (^p(Z))*, the 
dual space of ^p(Z) (in the sense e.g. of Kato [24], where the elements of the dual space are anti- 
hnear functional), it holds that A* : £9(Z) -^ £i{z) is the adjoint of A : iP{Z) ->■ iPiZ). Further 



\q- 



[24] A is invertible iff A* is invertible and, if they are both invertible, then \\A ^\\p = \\{A*) 
Similarly, A is Fredholm iff A* is Fredholm and, if they are both Fredholm then indA = —indA*. 

In this section we first compute the numerical range of the operator A'' in the case when b is 
pseudo-ergodic, which gives an upper bound on the spectrum S of A^. We then apply a variety of 
symmetry arguments to explore the relationship between spectral sets for matrices with one and 
two ±1 diagonals and between semi-infinite and bi-infinite matrices, and to explore the geometry of 
E and that of E^, the e-pseudospectrum of A^ on iP{Z) when b is pseudo-ergodic. Our final result 
shows that, roughly speaking, in the pseudo-ergodic case, the spectral sets are the same whether 
the matrix is semi-infinite or bi-infinite, and whether the matrix has one or two ±1 diagonals. 

These results are to some extent surprising: there is no expectation in general that the spectral 
sets associated with bi-infinite and corresponding semi-infinite matrices will be the same. A simple 
example is provided by the shift operator V-i. This is a Laurent operator (a bi-infinite Toeplitz 
matrix) whose spectrum is the unit circle and whose £'^ e-pseudospectrum is the e-neighbourhood of 
the unit circle. On the other hand the Toeplitz operator that is the shift operator restricted to £^(N) 
(a semi-infinite Toeplitz matrix) has spectrum that is the closed unit disc (e.g. [15]). An example 
closer to our case is studied in [39], where calculations are made of the spectra of random bi- 
diagonal bi-infinite and semi-infinite matrices, matrices which the authors term stochastic Laurent 
and Toeplitz operators, respectively, by which they mean a bi-infinite or semi-infinite matrix where 
each diagonal is either constant or has random entries, but with the random distribution constant 
along the diagonal. In the bi-diagonal case they study, which has the constant value 1 along the first 
superdiagonal and a random main diagonal, it is found [39] that the bi-infinite and semi-infinite 
matrices may or may not have the same spectra, this depending on the support of the probability 
density function for the random variables on the main diagonal. 

Our first result is a computation of the numerical range. By W{B) we denote the (2-norm) 
numerical range of the operator or matrix B, defined by (see Section 1.2) W{B) :— {(Bx,x) : 
\\x\\2 = 1}, where (•, •) denotes the standard i^ inner-product on C" or on £'^{S), with S" = Z or N, 
as appropriate. 

Lemma 3.1 For b e {±1}^, W{A''J C W^(A^) C W{A'') C A := {x + [y: x,y eR, \x\ + \y\ < 2}, 
and W{A^) ~ A if b is pseudo-ergodic. Similarly, W{A'i^) = A if b E {±1}^ is pseudo-ergodic, 

and E C A. 

Proof. For b £ {±1}^ and x G £'^{Z) with ||a;||2 = 1, defining a — xV-ix we see that 

{A''x,x) = ^{bk^iXk-i +Xk+i)xk = ^{bkO-k + ak) = ^[ak{l + bk) +i/3fe(l - bk)], 
fcez feez fcez 

where Uk = ^{ak) and Pk — ^(afe). Thus 

mA'>x,x)\ + \Q{A'>x,x)\ < ^{|afe|(l + 6,) + |/3fc|(l - bk)} < 2^ |afc| ^ 2||a||i. 

Now, since x S ^^(Z), x and V-ix must be linearly independent. Hence, by the Cauchy-Schwarz 
inequality, ||a||i = ||xF_ia;||i < ||x||2|lVlix||2 = 1. We have shown that W{A'') C A; it follows 
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that W{A'^J C VF(A^) C W{A'') from (19). From this it foUows, from standard properties 
of the numerical range (see the end of Section 1.2) that S C VF(A'') C A. But, since tti = 
[—2,2] U i[— 2,2] C S, this impHes that ±2 and ±2i are in W{A''), if b is pseudo-ergodic. Hence, 
if b is pseudo-ergodic, then, for every rj > there exist points r,s,t,u G W{A^) with |2 — r| < rj, 
|-2-s| <!], \2i~t\ <-q and |-2i-u| < 77. Since W{A^) is convex, this imphes that A C W{A^), 
and so W{A'') = A. A similar argument, using that E = specA^j. C VF(A'|_) if 6 G {±1}^* is 
pseudo-ergodic (that S = specA^ is established in Theorem 3.6 below), shows that M^(A'^) = A 
if 6 e {±1}^ is pseudo-ergodic. ■ 

Our next result elucidates the relationship between the spectral properties of matrices with one 
and two ±1 diagonals. One obvious symmetry result we use already in this lemma is that, since the 
coefficients 6, c e {±1}^ are real-valued, the spectrum and pscudospectrum of A^''^ are symmetric 
about the real axis. 

Lemma 3.2 For a,b,c G {±1}^, 

where d = aV-ia, so that 

spec A^-'^ = spec A^''^"'^ = spec A^^ spec,,, A^-'^ = spec,,, A'"^''''^ = spec,,, A^". 

Further, for A ^ spec,,, A''''^, ind (A'''^ — A/) = 0, and, for A ^ spccA'''^ and p E [l,oo], where 
q e [1, 00] is given by p^^ + q^^ = 1, 

WiA''- ~ A/)-i||p - IKA"'^-^'^ - A/)-i||p = ||(A^^ - A/)-i||p = ||(A^'^ - A/)-i||„ 

so that, for e > 0, 

spec^A^''^ = specPyl''''''^'' = spec^A''" = specf^'''". 

Moreover, for I < p < r < 2 and e > 0, 

spec^A^^= C specPA''-". 

Proof. For a,b,ce {±1}^, recalling (20) and noting that M^^ = Ma, 

In particular, choosing a so that d = c, this identity reduces to AfaA'''^Af^^ = A'"^, while, choosing 
a so that d = be, this identity reduces to MaA''''^M~^ — A'^'^ — {A^''^)* . The remaining results, 
except the last equation, follow since Ma is an isometric isomorphism, and using the properties of 
the adjoint listed immediately at the beginning of the section, and standard properties of Fredholm 
operators, e.g. [24, 23]. The last inclusion follows from the interpolation theorem of Riesz-Thorin, 
often called the Riesz convexity theorem [38, Chapter V, Theorem 1.3], which implies that, for 
A ^specA^''', 

\\{A'^' - \I)-% < maxiUA'-^- - A/)-i||„ [KA"- - A/)-i||,) = UA"-^- - A/)-i||,. 

■ 

Note that this lemma implies that, for l<p<2<(7<oo, where p"^ +9^^ = 1, 

spec^A^^^ C specPA'''^ = specfA'''^ C spec^A'''^ ^ spec^A^''^ 

In general, for a non-self-adjoint operator or matrix A, it need not hold that spec^A C spec^A for 
any distinct p,r £ [1, 00]. 

Exactly the same results hold in the semi-infinite case. Precisely, where M+ denotes the operator 
on £P(N) of multiplication by a e {±1}^, Lemma 3.2 holds also with Z replaced by N, Ma replaced 
by M^ , and all other operators replaced by their semi- infinite counterparts. Similarly, where Z?" is 
the diagonal matrix with the vector a — (ai, ..., a„) on the diagonal, the following finite dimensional 
version of the above lemma holds. 
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Lemma 3.3 For neN, ae {±1}", and b,ce {±1}""\ 

r^a Ab^c T^a Abd^cd 

^n^n ^n — ^n ; 

where d — {aia2, •■■, a„_ia„), so that 

specA^^'= = specAi^-'^'^ = spec A^^ 
Further, for p € [1, oo] and e > 0, where q G [1, oo] is given by p^^ + q^^ = 1. 

specPA'',f = specP^^'*'^'' = specP^,',^ = specf^^'^ 
Moreover, for I < p < r < 2 and e > 0, spec^A^'^ C specf Aj^'"^. 

A first application of the above lemmas is the following symmetry result (of. [22]). 

Lemma 3.4 For b G {±1}^, e > 0, and p G [l,oo], spec^'', spec^^gg^'', specf^'', specA^, and 
spec^A^ are invariant under reflection in the real and imaginary axes. Further, where S{b) denotes 
any one of these sets, it holds that S{—b) = iS{b). The set S, which is the set spec A*" — spec^gg A'' 
in the case that b is pseudo-ergodic, and, for e > and p G [l,cx)], the set T,p, which is the set 
spec I? A^ for b pseudo-ergodic, are invariant under reflection in either axis and under rotation by 
90°. 

Proof. We prove the results for A^ using Lemma 3.2; the proof for A^ using Lemma 3.3 is similar. 
That the entries of the matrix A'' are real implies the symmetry about the real axis. Defining 
a G {±1}^ by a^ = (—I)'', fc G Z, so that d = aV-ia is the constant sequence d = (..., —1, —1, ...), 
it follows from Lemma 3.2 that MaA^M~^ = —A^, which implies that the sets specA'', speCggg A*", 
and spec f A'' are also invariant under reflection in the origin, so that they are also invariant under 
reflection in the imaginary axis. Defining, instead, a G £°°(Z) by a^ = i*^, we obtain, similarly, 
that MaA''M~'^ = A'^''''^, where d = aV-ia so that dk = i. Thus MaA''M-^ = -i^-^ and we see 
that S{—b) ~ iS{b), where S{b) denotes one of spec A*", speC(,ssyl^, or specfA^. Where S{b) again 
denotes one of these sets, since b is pseudo-ergodic iff —b is pseudo-ergodic, that S{b) = S{~b) = 
iS{b) follows from Theorem 2.5. ■ 

The following lemma further elucidates the relationship between the spectral properties of semi- 
infinite and bi-infinite matrices. In this lemma for p E [1, oo] we let i^{Z) denote the closed subspace 
of odd elements of (P{2.), i.e. x G ^o{Z) iff x-k = —Xk, fc G Z, and let £f (Z) denote the closed 
subspace ofeven elements of^P(Z), i.e. a; G £f(Z)iffa;_fc = Xfe, fc G Z, so that F(Z) = £p(Z)®^p(Z). 
It is convenient to equip ^^(Z) with the norm ||a:|| := 2^^/P||a;||p, so that the extension operator 
E : £P(M) -^ iP{Z) given by {Ex)k = Xk, k € N, {Ex)o = 0, and iEx)_k = -Xk, fc G N, is an 
isometric isomorphism, as is the restriction operator P : £p(Z) — )■ £p{N) given by {Px)k = Xk, 
fc G N. (This change of norm does not effect the value of the induced norm of a bounded linear 
operator A on ^o(^)i ''^^'^ ^o does not affect the definition of specf A.) Further, let R : P'{'L) -^ £p{Z) 
be the reflection operator given by {Rx)k — x_i_fc, fc G Z. 

Lemma 3.5 Suppose b G {±1}^ with bk = I, k < 0, and let c = Kb. Then, for p G [l,oo], A^^'^ 
maps £P{Z) to l^CE) and maps ^f (Z) to £f (Z). Further, where A'^''^ denotes the restriction of A^''^ 

to ep{z), 

AY = EA\P. (29) 

Thus specA\ = specAj^-'^ C specA''''^ = specA'"^ and specfA+ = spedPAY C specfA^^'' = 
specfA'"^, for e > and p G [1, oo] . 

Proof. For x G £l{'L), {A^^''x)q = b^iX-i + CqXi = a;_i -h xi = and, for fc G N, (A'''=x)_fc = 
b-k-ix-k-i + C-kX-k+i = -CkXk+i - bk-iXk-i = -(A'''^a;)fc, so that A''^" : iP{Z) -)■ iP(Z). 
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Similarly, for x e ^e(Z) and fc e N, {J^''''x)-k = b-k-iX-k-i + c-kX-k+i = CkXk+i + bk-iXk-i = 
{A''-''x)k, so that A''-" : £^(Z) -^ £p{Z). Further, for k e N, {EA''_^x)k = fofc-iXfe_i + Xfc+i = 
bk-iXk-i + CkXk+i = {A\^'^x)k, so that (29) holds. Since E and P are isometric isomorphisms 
and E — P"^ , it follows that spec^^J. = spec^g''^ and that specfv4+ — specfAg'^, for e > and 
p € [1, oo]. The remaining results follow from Lemma 2.8 and Lemma 3.2. ■ 

Putting the results from the previous section and this section together gives the following 
characterisations of the spectrum, essential spectrum, and pseudospectrum in the pseudo-ergodic 



Theorem 3.6 // b,c,d E {±1}'*', e,f,g G {±1}^, and b, cd, e, and fg are pseudo-ergodic, 
then spec A^ = spec A^ = spec A'^ — spec A^'^ — spec^sg A\_ = speCg^g ^^ = spec^^g A*^ = 
speCgss A-'^^^ — E and, for e > and p e [l,oo], where q € [l,oo] is given by p^^ + q^^ — 1, 
specPA*; = specf A^'* = specf A'^ = spec^Af^s = SP = S?. Further, for I < p < r < 2 and e > 0, 

Y^l C Sf . 

Proof. From Lemma 3.2 and the remarks following that lemma we have that spec A'^^ — spec A'^f, 
speCgggA^ = speCgggA^, spcc^-'''^ = specA-''^, and speCggg A-^'^ = speCggg A-''^. From Theorems 
2.5, 2.9, and the remarks at the end of Section 2, we have moreover that if b and e are pseudo- 
ergodic then specA"^ = spec^gg ^"^ = speCggg yl'j_ = S. This implies that E C specA^, and that 
specA'j, C S follows from Lemma 3.5 which, together with Lemma 3.2, gives that specA!j_ C 
specA'''^'' — specvl^-'^'' C S. The results for the pseudospectrum are shown similarly, again using 
Theorems 2.5, 2.9, the remarks at the end of Section 2, and Lemma 3.2. ■ 



4 The relationship between the spectra and pseudospectra 
of finite and infinite matrices 

An obvious method to try to calculate the spectrum of an infinite matrix is to study the spectra 
of large finite submatrices of the infinite matrix and hope that these provide good approximations. 
In particular, one can apply this idea to the infinite matrix A^, and hope that the spectrum of 
the n X n matrix A^ , the intersection of the first n rows and columns of A^_^ , will approximate the 
spectrum of A^ well for n large. 

In general the spectrum and pseudospectrum of an infinite banded matrix may or may not be 
well-approximated by the spectra and pseudospectra of its finite submatrices (see [30] and the 
references therein for some discussion, with emphasis on the case of tridiagonal pseudo-ergodic 
matrices). In particular, there need be no relationship at all between the spectrum of a bi-infinite 
matrix and the spectra of its finite sections. A simple example is provided by the Laurent operator 
that is the shift operator V-i with matrix representation {aij)ij,£z, with Oij = Sij+i, whose 
spectrum is the unit circle. The Toeplitz matrices that are its n x n finite sections, (aij)i<i.j<„, 
clearly have zero as the only eigenvalue. 

The purpose of this section is to show that, for the particular class of pseudo-ergodic operators 
we are studying, there is a perhaps surprisingly close (given that our pseudo-ergodic operators are 
not self-adjoint or normal) connection between the spectral sets in the finite and infinite case. This 
connection is particularly close for the pseudospectra. 
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4.1 That the finite matrix spectral sets are contained in the infinite 
matrix counterparts 

For n e N, introduce the n x n matrices 

1 



In = 



and Jn — 



1 



so that /„ is the order n identity matrix. The proof of the foUowing result uses a similar construction 
to that of the bi-infinite matrix A^''^ in the proof of Lemma 3.5. 



Theorem 4.1 Ifb is pseudo-ergodic then, for n € N 
if- 



spec^°j C o-„ := [J specA^ C 7r2„+2 C spec A*" = E. 
/G{±i}"-i 





Figure 3: Our figure shows the sets an of all nxn matrix eigenvalues, as defined in Theorem 4.1, for n = 5, 10, ...,30. 
Note that in the first pictures (with only a few eigenvalues), we have used heavier pixels for the sake of visibility. 

Proof. If A e specA^, for some / g {±1}"^^, then Af^x — Xx for some non-zero x G C". Put 

X := JnX and An :— JnA(^Jn- Then 



AjiQ 



'Jn-^n'^n^n-^ 



'Jn-^n'^ 



Jri\X \X 
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and hence, using repeated reflections, i.e. by putting 

/ ■. 

1 



A 



c,d 







1 




1 


AL 


1 




-1 




1 



A^ 

-ii-71 



V 



/ : \ 



and X 



J 



V : / 



we get A'^^'^x — Xx with x E £°°{Z), so that A is an eigenvalue of A'^-'^ as an operator on ^°°(Z). 
Thus, applying Lemma 3.2 and Theorem 2.5, and noting that both c and d are periodic, with 
period 2n + 2, we see that A S spec A'^''' — spec A'^'' C 7r2„+2 C S = spec^''. ■ 




Figure 4: An illustration of the inclusion 0-4 C ttiq, which holds by Theorem 4.1. (The points in 0-4 are indicated 
by circled dots.) For similar figures for other values of n see [27]. 

In Figure 3 we plot the sets (t„, for n — 5, 10, ...,30 (note that each set cr„ is invariant under 
reflection in either axis or under rotation by 90°, by Lemma 3.4, and see [5] for smaller plots of 
these sets for n = 1, ..., 30). By the above theorem, (t„ C Tr2n+2 for each n, so that 



CToo := [J CTn C TTo 
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The inclusion (t„ C 7r2„+2 is illustrated for n = 4 in Figure 4. 

An interesting question, alluded to already in Section 2, is whether tToo, which is contained in S, 
or CToo, which is a countable subset of tToo, are dense in S, the spectrum of A^ for b pseudo-ergodic. 
Of course, we do not know what S is, so that this question is difficult to resolve! We do know 
however (Theorem 2.7) that the unit disc D C S, and we can consider the question as to whether 
TToo or tJoo are dense in D. Recall that the sets 7r„, for n = 5, 10, ..., 30, are plotted already in Figure 
2. Studying Figures 2 and 3, it appears that there is a "hole" in both cr„ and 7r„ around the origin, 
though these holes appear to be reducing in size as n increases. And in fact, as mentioned already 
in Section 2, it has been shown recently that tToo is dense in D. Further, it appears to us plausible, 
comparing the two figures, to conjecture that CToo is dense in tTqc and so dense in D. 




■ U 






' » 




. .-■«'■■.: 



■(lb * 



Figure 5: This is a zoom into 0-25 - the 5th picture of Figure 3. The location of this zoom is near the point 1 + i, 
which is the midpoint of the northeast edge of the square W{A^) = A. The picture clearly suggests self-similar 
features of the set (T25 ■ 

Figure 5, taken from [5], zooms into the part of the set (T25 around 1 + i. Intriguingly this set, 
the collection of all eigenvalues of a set of 2^^ matrices of size 25 x 25 (25 x 2^* = 419, 430, 400 
eigenvalues in all!), appears to have a self-similar structure. We have no explanation for these 
beautiful geometrical patterns, and it is not clear to us how to gain insight into the geometry of 
this set. 

In the next theorem and corollary we show the analogue of Theorem 4.1 for pseudospectra. 
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Theorem 4.2 // b is pseudo-ergodic and n G N then, for all X E C \Y,, f E {±1}" ^, and 
pe [1,00], 

\\{Ai-XIr.)-% < \\{A'-XI)-%. 

Proof. Let A G C \ S and / e {±1}"^\ so that A^ - XIn is invertible by Theorem 4.1, and let 
p e [1,00]. Put M := \\{Al - A/„)^^||p. For every 6 > 0, there exists an cc = (a:i, ...,x„)^ S C" 
such that ||a;|| = 1 and y := {Af^ — A/„)a; has ||y|| < jfzrg- Now let c,d & {±1}^ be the sequences 
in the matrix A'^''^ introduced in the proof of Theorem 4.1. At this point our current proof has to 
bifurcate depending on the value of p. 

Case 1 : p = 00 

Define x G £°°{Z) exactly as in the proof of Theorem 4.1. Then y := {A'^-'^ — XI)x is of the form 

y = (• • • ,2/T, 0, (J„2/)T, 0,2/T, 0, (Jny^ , • • 0^ e ^°°(Z) and ||y||^ = ||y||^, as well as \\x\\^ = ||x||^, 

so that 

IItII IItI 

" "°°"lly|loo bIL 

Case 2 : p < 00 

For any m E N, let ■r'^™^ be the sequence x from case 1, but with all entries of index outside 
{— m(n + 1), . . . , m{n + 1)} put to zero (where we suppose that the sequence x is numbered so 
that at index zero there is one of the entries between x and x of x, so that Xq = 0). Then 
~(m) ._ ^^c,d _ xj'j'^im) jg ^Yie same as y from Case 1 for entries with index between —m{n+ 1) + 1 

and m{n + 1) — 1, is zero outside {—m{n + 1) — 1, ... , m{n + 1) + 1} and we have y_^/„ , j^-, = xi 

if m is even and y , , -.^ — Xn if to is odd, while ir, , ,-, = —x^ if ra is even and y , , ■,-, — —Xr, 
if m is odd. As a result, we find that 



p('")||^ = 2TO|la;||^ and Hy^™)]!^ = 2to ||y||^ + | 



2 |xi|^ if m is even, 
2 \x„\P if m is odd. 



From ||a;||p = 1, ||y||p < jj^ and |a;i|,|a;„| < ||a;||p = 1 we hence get that Hx^^^Hp — 2m and 
||y(™)||P<2TOp7^+2, sothat 



erf ^n-l„P ^ ll^'™'llp . 2to 



wiA^-'^ - xirx > ^(^ > 

\\y^ '\\p 



(M-Sy ' ^ (M-S)v ' m 



In either case, Case 1 or 2, these inequalities hold for all (5 > and all m € N. Hence, and 
applying Lemma 3.2, 

\\{Ai - xi^rx = M < \\{A^''' xirx = \\{A^' ^irX < ii(^' - ^irX^ 

where the last inequality follows by [8, Theorem 5.12(ix)], since A'^'^ is a limit operator of A^ by 
Lemma 2.4. ■ 

The following corollary is immediate. 
Corollary 4.3 // b is pseudo-ergodic and n E N then, for all e > and p G [1, 00], 



7P •= 



U spec^A^ C specP^'' = SP, 

CG{±1}"-1 



and in particular specP A'^ C spec^A''. 
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4.2 Convergence of the finite matrix spectral sets to their infinite matrix 
counterparts 

As we have remarked at the beginning of this section, it is not clear that the spectrum of a 
general banded matrix should have anything to do with the spectra of its finite submatrices. In 
particular, it need not be the case either that the spectrum of a large finite submatrix is contained 
in a neighbourhood of the spectrum of the corresponding infinite matrix, or that the converse 
statement is true. But the situation is somewhat more positive for the pseudospectrum, namely 
that, as we show for a general tridiagonal matrix as our first result of this section (and our method 
of argument applies to banded matrices more generally), the e-pseudospectrum of the infinite 
matrix is contained in the e'-pseudospectrum of an appropriately chosen n x n submatrix, for a 
given e' > e, provided n is sufficiently large. The argument is based on a standard and rather 
obvious idea: the point is that every eigenvector, or approximate eigenvector, of the infinite matrix 
is, when truncated in a careful way, also an approximate eigenvector of the finite matrix. 

The opposite statement is, in general, false; an approximate eigenvector of a large finite matrix 
is an approximate eigenvector also of an infinite matrix, but the infinite matrix that it is an approx- 
imate eigenvector of need not be the infinite matrix whose spectrum one wishes to approximate! 
(One recent result which expresses this idea very precisely in the l'^ case for a version of the finite 
section method for the class of general pseudo-ergodic tridiagonal matrices is [30, Theorem 2.14].) 
But, for the pseudo-ergodic operators A^ and A^ that we are studying, we have also shown, in 
Corollary 4.3 and Theorem 3.6, that spec^ A'^ C spec^A^ — spec^A^. Putting this result together 
with Theorem 4.4 proves that spec^A^ /• specfA'' (using the notation of the introduction and 
Section 1.2). 

Theorem 4.4 Suppose that A = (aij)j jez is a bi-infinite tridiagonal matrix with M :— supj \aij \ < 
CXI. Define the semi-infinite matrix 74+ by A^ = (aij)i.jgN and, for £,m Cz N with i < m, define 
the finite matrix of order m + 1 ^ i by Ai^m = io-ij)i,je{i....,m}- Then, for every e' > e > and 
p G [1, oo], there exists N gN such that 

spec^A C spec^, Ag^rn, for £ < -N and m> N, (30) 

and 

spec^^+ C spec^,Ai,,„, for m> N. (31) 

Proof. We will prove (30). The proof of (31) is similar. 

As a first step we will show that, given some £ > and p G [1, oo], for every A S spec^A there 
exists N G N (depending on A) such that A G specf A^^m if ^ < —N and m > N. We then combine 
this result with a compactness argument to obtain the proof of the theorem. 

So suppose that e > 0, p G [l,oo], and that A G specf^. Then, by Theorem l.l(iv), either 
A G specP„i„j_^A or A G spec^^j^^ ,^A*, where p-^ + q-^ = 1. 

Suppose first that A G s>pec^^^^^ ^A, i.e. that there exists x G iP{Z) with H^llp — 1 and e := 
\\y\\p < e, where y :— {A ~ XI)x. In the case p < cx), let x := {xi, ...,Xm)'^ and y := (^^,m — 
XIm+i-e)x, so that yk — yk, k — £ + 1, ...,m — 1. Since \xk\ — )■ as |fc| — >• oo, it is easy to see 
that we can, given 6 > 0, choose N such that \\y\\p < e + 5 and ||a;||p > I — S whenever £ < —N 
and m > N. But this implies that A G spec^A^^m if N is large enough and £ < —N and m > N. 
In the case p = oo we have to modify this argument slightly. Given £ < —N and m > N put 
X — (xe, ...Xm)'^ with Xk ■= ujN{k)xk, k — £, ..., m, and WAr(fc) -.^ max(0, 1 — |fc|/(A^ — 1)), fc G Z, 
and let y :— (A^.m ^ ^Im+i-i)x. Then, for i — £, ..., m, 

\yi\ = |aij-i^j-i + a.jXj + a,jj+i-Ej+i| 

= \ujN{j)yj + aij^i{u>NiJ - 1) - '^nU))^]-! + ai,j+iiujN{j + 1) - ujn{J))xj+i\ 
< \y,\ + 2M\\x\U/{N - 1) < £ + 2M/{N - 1), 

23 



since \ujn{J) — ujn{J + 1)| < (^ — 1) "^> for k £ Z. Since also, for each k £ i, ...,m, xj. — > Xk as 
iV — > oo, it is clear that, for every 5 > 0, ii N is chosen large enough, then ||x||oo > ^~ S, and also 
||y||oo < £ + S. But this implies that A £ specf ^f,m if N is large enough and i < —N and m> N. 

If A e speCpQjjj^ ^A* then essentially the identical argument shows that A e spec|A|„j. But this 
implies that A G specf A^.m [40, Section 4]. This completes the proof of the first step. 

To finish the proof of the theorem we argue as follows. Given e' > e > and p £ [l,oo], let 
T] := {e'~e)/2, and e* = e + 77. Let S := specf^, and let O := {A + 77D : Xe S}. Then O is an open 
cover of the compact set S, and so has a finite subcover, i.e. there exists a finite set A C specfA 
with 5 C UAeA('^ + ^''^) ~ rjD + A. Now A C specf A C spec^. A. Applying the result shown in the 
first step, we see that we can choose N so that, for £ < —N and m > N, A C spec^.A^^m. Thus 
specPA cScr]D + Acr]B + spec^. A^,™ c spec^,A^_„, by (15). ■ 

To apply this result, for £,m £ Z with £ < m, let A'^ ^ denote Ai,rn, the matrix of order m + l — £ 
as defined in the above theorem, in the case that A — A^. So, in particular, A\ „ = Af^ for n E N. 

Corollary 4.5 If b E {±1}'^ is pseudo-ergodic then, for every e > and p £ [l,cx)], 

spec^A^ = specPAi_„ /^ spec^A''^ = E^, as n ^ 00. 
If b E {±1} is pseudo-ergodic then, for every £ > and p G [l,oo], 



specPA^ „ /" specPA'' = E^, as £ ^ - 



00 ana to — > 00. 



Proof. We will prove the second of these statements. The proof of the first is similar. From 
Corollary 4.3 and Theorem 4.4, given any e' G (0, e) there exists N E N such that 

spec^,^ C spec^ Aim C spec^v4, for £ < -N and m> N. (32) 

Since, from (18), spec^,yl /* spec^A as e' — > e~, it follows that specf A^ ^ /* specfA'' (which is 
equal to E^ by Theorem 3.6), as ^ — > —00 and tti — > 00. ■ 

A similar result holds for the convergence of the numerical range, as an instance of the general 
result Theorem 1.3. Note that while convergence of the pseudospectra needs that b is pseudo- 
ergodic, to ensure that the matrix pseudospectra are contained in the operator pseudospectra, the 
corresponding inclusion (19) for numerical ranges holds for any bounded linear operator, so that 
we need no constraint on b. The following is thus an immediate corollary of Theorem 1.3 and 
Lemma 3.1. 

Corollary 4.6 If b e {±iy'^ then 

WiAi) = WiA\j / W{AX), as n ^ 00, 

with W{A'^) = A ifbis pseudo-ergodic. If b £ {±1}^ then 

W{A\„;) /"WiA^), as £ ^ -00 and m ^ 00. 

with W{A'') = A if b is pseudo-ergodic. 

4.3 Quantitative convergent approximations to the spectrum and pseu- 
dospectrum 

In this section we present numerical algorithms for approximating S and S^ which are, respectively, 
from Theorem 3.6, the spectrum and the £^ e-pseusdospectrum of both A'' and A'j_ in the case 
when 6 is pseudo-ergodic. 
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The previous subsection already provides potential methods for computing these sets. We have 
that, if 6 e {±1}" is pseudo-ergodic, then 



E^ = lim spec^^^ (33) 



This then implies, by (18), that 



S = lim lim spec?A^. (34) 

In principle, these equations can be used as the basis of algorithms for computing S^ and S. In 
particular, to approximate S^ one uses the sequence of sets spec^A^, n — 1,2, ..., which can be 
computed as described in Section 1.2. The difficulty with this scheme is that one has no idea of the 
rate of convergence of spec^A'^ to E. Indeed it is clear that it can be arbitrarily slow: to see this 
consider that if c e {±1}'*' is pseudo-ergodic, then so is 6 e {±1}^ if &„ = c,„ for all sufRciently 
large m. But this means that it can hold that b is pseudo-ergodic and that bm = 1, for 1 < m < N, 
with N arbitrarily large. If this is the case then A^^ is self-adjoint and thus, and by Lemma 3.1, 
specA^ C (-2,2) and spec^A^ = specAf, + eD C (-2,2) + eD, for n < N. So if, e.g., N = 10^ 
then, while ultimately spec^A^ -^ S^, there is no early sign of this. 

The situation with (34) is rather worse. This equation implies that there exists some sequence 
of positive reals e„ — )■ for which it holds that 

spec^^A^ -> E, 

but provides neither a recipe for choosing the e„ nor any guarantee of the rate of convergence. 

The source of the difficulty regarding the rate of convergence can be traced back to Theorem 
4.4 and its proof, this theorem a key ingredient in the proof of Corollary 4.5 and so of (33). This 
theorem guarantees that, for every e' > 0, S^ = spec^A^ C spec^iA^ for all n sufficiently large, 
but gives no idea of how large n should be. And indeed we have argued above that there is no 
upper bound on how large n may need to be for this equation to hold for a given pseudo-ergodic b. 

This difficulty has been resolved in recent work by the authors [6], who quantify, for general 
tridiagonal matrices, by a sharpened version of the arguments of Theorem 4.4, adapted particularly 
to the case p = 2, exactly how e' should depend on n in (31), but at the expense of replacing in this 
equation the pseudospectrum of a single n x n submatrix by the union of the pseudospectra of all 
possible nx n principal submatrices. The results in [6] are much more general, but we will restrict 
the exposition here to how these results apply to the bi-infinite matrix A'' with b G {±1}'^. Using 
the notation of Corollary 4.5, the result shown in [6] (or see [10, Corollary 3.7]) is the following 
when applied to A'^: 

Theorem 4.7 For b G {±1}^, £ > 0, and n e N, 

spec^A'' c IJ spec2^^^A^^^+„_i, 

where Sn = 4sin0„ < 2n/{n + 2), with On the unique solution in the interval —. — -, — -- 

of the equation 2 cos ((n + 1)0) ~ cos ((f^ — 1)^) . Further, 



spec A'' c U specie Al^^^_^ = |J Spec^^vl'' 



e.,i+n-i- 
1^1 



An important point is that the unions of pseudospectra over i gTL'vh the above equations reduce 
to finite unions, because there are only 2"^^ distinct nx n matrices A^ with c e {±1}"^^. In 
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the notation introduced in Corollary 4.3, it must hold that [J^^^speciA'^ ^^^_-^ C c^^, for every 
77 > 0. For small values of n, e„ in the above theorem can be calculated explicitly, in particular 

ei == 2 and £2 = V2. (35) 

Example 4.8 As a first example of application of the above theorem, consider the case when 
bjn ^ I for each m. Then A^^_|_„_j^ = Aj „ = A^ for each i. Further, this matrix is self-adjoint, so 
that spec^A^ = specA^^ + ryD, for every 77 > 0. Thus the statements of the theorem reduce to 

spec A^ C spec A^ + e„D and spec^A^ c specAf^ + (e + e„)ID>, e > 0. (36) 

In this simple case we can compute the above sets explicitly, to check that the above inclusions hold, 
finding that specA'' = [-2,2], spec^A'' = [-2,2] + eD, and specAf^ == |2cos;^ : j ^ l,...,n|. 

Elementary calculations show that the inclusions (36) do hold in this case, in fact one can calculate 
(see [10, Section 3.2.2] for details), if £„ were replaced with ej^ < e„ in the above inclusions, the 
smallest value of e*^ for which the inclusions would still hold. This is e^ = 2, e* = 2 sin(7r/(2(n+l))) 
if n is even (in particular £3 = l)i ^nd e*j = sin(7r/(n + 1)) if n > 3 is odd. Thus Sn/^n ~ 1 '^or 
n — 1 (the bound (36) is sharp for n = 1) and en/s'n —^ 2 a.s n —> cx). D 

The main example of interest to us here is the case where b e {±1}^ is pseudo-ergodic. Recall 
from Theorem 3.6 that spec^'' — S and spec^A'' = E^ in that case. Combining Theorem 4.7 with 
Theorem 4.1, Corollary 4.3, (13), (18) and (17), we obtain the following result. 

Theorem 4.9 For e > and n £N, 

cr„ C S C 0-2 ,^^ C Sl^ and E^ C cr,le+£„ C ^l+e„ , (37) 

where £„ is defined as in Theorem 4-7. Further, 

CT„ + eD C E^, for £> 0, and E^_g^ C cr^^^ C E^, for e > En- (38) 



Moreover, as n -^ 00, a^ ^ \^ and, for e > 0, a'^^ ^^^ \ E^ and cr^ ^ /* E^. 

In most respects this result is superior to Corollary 4.5. It provides both upper and lower 
bounds for E^ , moreover these converge to E^ as n — >■ 00 at guaranteed convergence rates (at least 
as fast as Eg-^^ and E^+g^^, respectively). Further, the theorem provides an upper bound which 
is convergent to E, at least as fast as E^ . Of course, that the convergence rates are guaranteed is 
at a cost: evaluating a'^ for some n G N and 77 > requires exponentially large computation for 
n large. Precisely, using the characterisation (11), we see that 

al,, = {A e C : 5„(A) < 77}, (39) 

where 

SniX):= min s^i„(A^ - A/„), A G C. (40) 

ce{±i}"-i 

Clearly, computing S'„(A) for a particular A, to check membership of cr^ requires calculation of 
the smallest singular value of 2"~^ matrices of order n. Note that it follows from (12) that 

|^„(A)-5„(M)|<|A-/i|, A,MeC. (41) 
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Figure 6: Plots, for n = 6, 12 and 18, of the sets a'^ ^ , which are inclusion sets for S = spec A'', when b S {±1}^ 
is pseudo-orgodic. Also shown, overlaid in red, is the square A, with corners at ±2 and ±21, which is W{A''), the 
numerical range of A . Overlaid on top of that in blue is the set ttso U B which, by definition and Theorem 2.7, is a 
subset of S. 

In Figure 6 we plot a'^ ^^, for n = 6, 12, and 18. Each of these sets contains E, by Theorem 4.9, 
and note that each set is invariant under reflection in either axis or under rotation by 90", by Lemma 
3.4. On the same figure we plot the square A which, by Lemma 3.1, also contains S. It appears 
that, for n < 18, A C cr^ J .If this were to hold for all n £ N then it would follow, from Theorem 
4.9, which tells us that u^ ^ \ E, and Lemma 3.1, which tells us that E C A, that E = A. It 
seems impossible from these plots to take an educated guess as to whether or not A C cr^ ^ holds 
for all n, not least because the convergence rate of a^ ^ to E may be slow: Theorem 4.9 tells us 
that dist(CT2^, E) < dist(E2_^, E) but it follows from (13) that dist(E2^, E) > £„ « 27r/(n + 2). 

We have not been able to produce similar plots to those in Figure 6 for much larger values of 
n because of the large computational cost. But it is feasible to compute iS'„(A) for a single A for 
larger n. We have carried out this computation for A = 1.5 + 0.5i, a quarter of the way along one 
of the sides of A. Computing in standard double-precision floating point arithmetic we find that 

S'34(l-5 + 0.51) = 0.17201954132506... > £34 = 0.169830415547956... . (42) 

This implies that 1.5 + 0.51 ^ ct|4 ^ ^ and so 1.5 + 0.51 ^ E, which of course implies that E is a strict 
subset of A. In fact, in view of (41) and the symmetries of E noted in Lemma 3.2, the inequality 
(42) implies more, namely that 

(±(1.5±O.5i) + 77D)nE = 0, 

for 77 == £34 - 534(1.5 + 0.51) = 0.0021891257771.... 

We note that the computation required to evaluate £'34(1.5 + 0.51) and so establish that 1.5 + 
0.51 ^ E is considerable: we need to evaluate the smallest singular value of 2'^'^ « 8.6 x 10^ matrices 
of order 34 (of course these computations are ideally suited for parallel implementation). We 
note that it seems to be necessary to use n as large as 34, in that other computations show that 
533(1.5 + 0.51) < £33, so that 1.5 + 0.51 e crfg^^^^. 

5 The Random Case and Concluding Remarks 



We finish this paper by spelling out the implications of the above results for the finite matrices A. 

tb.c 



b 



and A^'^, the bi-infinite matrices A^ and A'^-'^, and the semi-infinite matrices A'', and A_i_'^, in the 



case when the entries of b and c are random variables taking the values ±1. 
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Theorem 5.1 Suppose that the entries of b £ {±1}^ are iid random variables, with Pr(6m = 1) € 
(0,1). Then: 

(i) specA^ c S, specA*^ c S, with speCgssA^ — specA'' = speCggg^!^ = specAj^ = S almost 
surely. 

(ii) T4^(A';_) C VK(A'') C A, with W{A^) = W{A\) = A almost surely. 

(iii) For n e N, spec Aj^ C S and VF(A^) C A, and, as n -> oo, VK(^fJ /^ A, almost surely. 

(iv) For e > and p e [1,00], specfA*" C Sf, specfyl'; C Sf, wii/i specfA'' = specfA^ = S§ 
almost surely. 

(v) For £ > 0, p e [1,00], and n G N, specfA^ C S^ and, as n ^ 00, specfA^ /* Ef, almost 
surely. 

Similarly, if b,c € {±1}^, and t/ie entries of be are iid random variables, with Pr(6„iC„i = 1) S 
(0, 1), then (i)-(v) hold with A^ , A^j^, A^, replaced by A^''^ , A_^ , and A^f- , respectively. 

Proof. To see that (i)-(v) hold, note that, by Lemma 2.3 and the remarks at the end ol Section 
2, the condition of the theorem imply that b and also 6^ := (&i,62i---) are pseudo-ergodic with 
probability one. Then (i) follows from the definition of E in Theorem 2.5, and from Lemma 3.5 
and Theorem 3.6. That (ii) and (iii) hold follows from Lemma 3.1, Theorem 4.1 and Corollary 
4.6. That (iv) holds follows from the definition of Ef in Theorem 2.5, and from Lemma 3.5 and 
Theorem 3.6. Finally, (v) follows from corollaries 4.3 and 4.5. That (i)-(v) hold for the case where 
A**, A^, A^ are replaced by A**'^, A^^, and Aj^'^, respectively, and the entries of be are iid random 
variables, with Pr(6m Cm = 1) € (0, 1), follow using the same results, with the help of lemmas 3.2 
and 3.3 and the observations between these lemmas on the semi-infinite case, and noting first that 
these assumptions imply that be and (&iCi, 62C2, •.•) are pseudo-ergodic, with probability one. ■ 

Of course, in the above theorem A = {z = a + i6 : a, 6 e M, |a| + |6| < 2} and E and E^ are 
as defined in Theorem 2.5. The following theorem summarises, for the convenience of the reader, 
what we have established in the sections above about the compact set E and the bounded open sets 
E§. Recall that CToo and tToo are defined by (7) and (6), respectively, cr^ ^ is defined in Corollary 
4.3 (and see (39) for p — 2), e„ is defined in Theorem 4.7, and S'„(A) in (40). 

Theorem 5.2 For e > and p £ [1, cx)], where q G [1, 00] is given by p^^ + q^^ = 1; 

(i) D C E C A, and E is a strict subset of A provided £34 — 6*34(1.5 + 0.5i) > 0, for which see 

(vi). 

(ii) tToo C TToo C E C Ef C E^, , for e' > e. (See Figures 2, 3, 4 and 5 for visualisations of an and 
7r„, for n G N, and their interrelation. ) 

(iii) E and E^ are invariant under reflection in the real and imaginary axes and under rotation 
by 90°. 

(iv) Ef = E| and E§ C S- z/ 1 < r < p < 2, so that E? = n.^n,,,, EJ. 



(v) As n ^ 00, CT^ ^ \ E, <Jn e+e \ ^e' '^'^'^ '^ne /* ^e ' (^^^ Figure 6 for visualisations of 
al,^,forn^'6" 12,18.) 

(vi) For X ~ 1.5 + 0.5i, provided rj = £34 — S'34(A) > (and floating point calculations give 
ri w 0.00219;, it holds that ±(1.5 ± 0.5i) + r;D n E = 0. 
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Proof. Part (i) follows from Theorem 2.7 (taken from [5]) and Lemma 3.1, and that E is a strict 
subset of A holds, as discussed at the end of 4.3, provided 77 = £34 — 5*34(1.5 + 0.51) > 0. Part (ii) 
is Theorem 4.1, with Sf C S^, because Sf = spec^A'' if 6 e {±1}^ is pseudo-ergodic (Theorem 
3.6). Part (iii) is Lemma 3.4, (iv) is from 3.2, (v) is part of Theorem 4.9, and (vi) is from the end 
of Section 4.3. ■ 

It is clear from the above results that we understand well, in Theorem 5.1, the interrelation 
between the numerical ranges and pseudospectra of the semi-infinite, bi-infinite, and finite random 
matrix cases, and have shown that the almost sure spectrum is the same set S for the semi-infinite 
and bi-infinite cases, and contains the spectrum in the finite matrix case. Interesting open questions 
are whether or not, similarly to the analogous results for the pseudospectra, specA'^ /^ Y, almost 
surely as n — >■ 00, which would imply that cr 00 is dense in S, so that tTqc is dense in E. (That (Too 
is dense in S was conjectured in [5].) Note that, if it does hold that spec Aj'^ /* S almost surely, 
then both Figs 2 and 3 are visualisations of sequences of sets converging to S. 

Regarding the geometry of S (and of the pseudospectra S^), we have some information in 
Theorem 5.2, including in the last part of this theorem establishing a computable sequence of sets 
converging from above to S (a sequence of three of these plotted in Figure 6). However there is 
much that is not known. Is E connected (which would imply, by general results on pseudospectra 
[40, Theorem 4.3], that also Ef is connected)? In fact, is E simply-connected? What is the 
geometry of the boundary of E, and the geometry of the sets (T„, the finite-dimensional analogues 
of E (cf. Figure 5)? We have conjectured in [5] that E is a simply-connected set which is the closure 
of its interior and which has a fractal boundary, which is plausible from, or at least consistent with. 
Figure 6, if it holds that W^ — E. Our methods and results provide no information about what is 
a usual concern of research on random matrices, to obtain asymptotically in the limit as n — ^ 00 
the pdf of the density of eigenvalues, except, of course, that we have shown in Theorem 5.1(iii) 
that the support of this pdf is a subset of E. 

There are many possibilities for applying the methods introduced in this paper to much larger 
classes of random (or pseudo-ergodic) operators. For some steps in this direction we refer the 
reader to [30, 6, 9]. 
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